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Abstract 

We present a preliminary measurement of the time-dependent CP asymmetries in decays of 
neutral B mesons to the final states -D* =F 7r ± , using approximately 82 million BB events collected by 
the BABAR experiment at the PEP-II storage ring. Events containing these decays are selected with 
a partial reconstruction technique, in which only the high momentum and the low momentum 
pion from the D*^ decay are reconstructed. The flavor of the other B meson in the event is 
tagged using the information from kaon and lepton candidates. We measure the time-dependent 
CP asymmetry A = —0.063 ± 0.024 (stat.) ± 0.017 (syst.). We interpret these results in terms of 
the angles of the unitarity triangle to set a bound on | sin(2/3 + 7)|. 
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1 INTRODUCTION 



Measuring the angles of the Cabibbo-Kobayashi-Maskawa (CKM) unitarity triangle [1] will allow 
us to overconstrain this triangle and to test the Standard Model interpretation of CP violation 
in the quark sector. A crucial step in this scientific program is the measurement of the angle 
7 = arg {-V ud V* b /V cd V* b ) of the unitarity triangle. 

The neutral B meson decay modes B° — > D*^^, where h is a light hadron (ir,p, a±), have 
been proposed for use in measurements of sin(2/3 + 7) [2], where [3 = arg (—V c dV* b /VtdV t l) . Since 
the time-dependent CP asymmetries in these modes are expected to be of order 2%, large data 
samples and multiple decay channels are required for a statistically significant measurement. The 
technique of partial reconstruction of D*~ mesons, in which only the soft (low momentum) pion tt s 
from the decays D*~ — > D°ir~ or D* + — ► D°7r+ is reconstructed, has already been used to select 
large samples of B meson candidates [3] . 

This paper reports the preliminary results of a study of CP-violating asymmetries in B° — > 
Z?* =F 7r :l= decays using the partial reconstruction technique. The analysis procedures for the selection 
of the signal and the reconstruction of the decay-time difference between the two B mesons in the 
event are essentially the same as those we have already applied to the measurement of the B° 
lifetime [4]. 



2 PRINCIPLE OF THE MEASUREMENT 

The decays B° -► D* T tt^~ may proceed via a d -> cud or a b — > ucd amplitude (see Figs. 1 
and 2). Interference between these amplitudes through B° — B° mixing produces time-dependent 
CP-violation observables [2] . The probability that a state produced at time as a B° or B° decays 
into the final state D* T 7r ± at time t is 

P(B° -► D^TT^t) = -^ e ~ Wr [1 ± C cos(Am d t) + S T sin(Am d t)] , (1) 

P(B° -» D* T 7r ± )(t) = ^~^ T [1 =F Ccos(Am^) - S T sin(Am d t)] , (2) 
where r is the B° lifetime, Am^ is the B° — B° mixing frequency, and we have defined 

C = 



1-r 2 

1+r 2 ' 



S ± = Y ^_ s in{2f3 + 1 ±5). (3) 

Here 5 is an unknown CP-conserving phase 7 , and r is the ratio of the magnitudes of the b — > ucd 
and b — ► cud amplitudes 

r \A(B^D*-n+)\ 

\A(B° D*+tt-)\' K ' 

Since the b — > ucd amplitude is doubly Cabibbo-suppressed with respect to the 6 — > cud amplitude, 
one expects r w 2%. Due to the small value of r, we use the approximations 

C « 1, 

5 ± w 2rsin(2/3 + 7 ±(5). (5) 



7 The definition of 5 is subject to additional 7r terms [5] that we ignore, as they are redundant with the discrete 
ambiguity 2f5 + 7 — > 2/3 + 7 + 7r, <5 — > 5 + 7r. 
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In principle, r can be measured from the first two terms in the square brackets of Eqs. 1 and 2. In 
practice, this requires sensitivity to terms of order r 2 , which is not available with the statistics of 
current data sets. 




Figure 1: Feynman diagram for the Cabibbo- Figure 2: Feynman diagram for the Cabibbo- 
favored decay B° — ► D*~it + . suppressed decay B° — ► D*~ir + . 

Due to the small value of r, a large number of signal events is required in order to observe 
and measure the small time-dependent CP asymmetry. In the partial reconstruction method, the 
decay B° — ► D* t tt^ is identified by reconstructing only the hard (high momentum) -k^ and the 
soft pion irj from the decay of the D*^. The four-momentum of the unreconstructed neutral D 
meson produced in the D*^ decay is calculated from the two observed tracks and the kinematic 
constraints relevant for signal decays. Partial reconstruction provides a way to obtain very large 
single-mode signal samples, by making use of events that cannot be fully reconstructed. 

3 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring. 
The data sample consists of 76.4 fb _1 collected on the T(4S) resonance, and 7.6 fb _1 collected at 
an e + e~ center-of-mass (CM) energy approximately 40 MeV below the resonance peak. Samples 
of simulated events with an equivalent luminosity four times larger than the data were analyzed 
through the same analysis chain. 

The BABAR detector is described in detail elsewhere [6]. We provide a brief description of the 
main components and their use in this analysis. Charged particle trajectories are measured by a 
combination of a five-layer silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) in 
a 1.5-T solenoidal magnetic field. Tracks with low transverse momentum can be reconstructed 
in the SVT alone, thus extending the charged-particle detection down to transverse momenta of 
~ 50MeV/c. Photons and electrons are detected in a CsI(Tl) electromagnetic calorimeter (EMC), 
with photon energy resolution a E /E = 0.023(E/ GeV)" 1 / 4 0.019. A ring-imaging Cherenkov 
detector (DIRC) is used for charged particle identification, augmented by energy loss information 
from the SVT and DCH. The instrumented flux return (IFR) is equipped with resistive plate 
chambers to identify muons. 

4 ANALYSIS METHOD 

4.1 Partial Reconstruction of B° — > £)* =F 7r ± 

In the partial reconstruction of a B° — > D*^^ candidate, only the hard pion track from the B 
decay and the soft pion track tt s from the decays D*~ — > D°irJ or D* + — ► D°irf are reconstructed. 
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The cosine of the angle between the momenta of the B and the hard pion in the CM frame is then 
computed: 

M 2 D ^-M 2 B0 -M^ + E CM E h 
cos 6 Bh = — ^ — |— , 6) 

where M x is the nominal mass of particle x [7], Eh and ph are the measured CM energy and 



momentum of the hard pion, Eqm is the total CM energy of the beams, and ps = y E^ M /4 — M^ . 
Events are required to be in the physical region \ cos 9 Bh\ < 1. Given cos 8 bh and the measured 
momenta of the 7T/j and ir s , the B four- momentum can be calculated up to an unknown azimuthal 
angle <\> around ph- For every value of 0, the expected D four-momentum Vd(<P) is determined 
from four-momentum conservation, and the corresponding ^-dependent invariant mass m((f>) = 
\/\Vd{4>)\ 2 is calculated. We define the missing mass m m i ss = ^ [m max + m m i n ], where m max and 
w-min are the maximum and minimum values that m(0) may obtain. In signal events, m m i ss 
peaks at the nominal D° mass M D o, with a spread of about 3 MeV/c 2 . The distribution for 
combinatoric background events is significantly broader, making the missing mass the primary 
variable for distinguishing signal from background, as described below. We use four-momentum 
conservation to calculate the CM D momentum vector with the arbitrary choice = 0, and use 
this variable as described below. 



4.2 Backgrounds 

In addition to B° — ► D*^tt^ events, the above procedure yields a sample containing the following 
kinds of events: B — ► D* T p =t ; peaking BB background (other than B — ► D*^p^), defined as 
pairs of tracks coming from the same B meson, with the soft pion originating from the decay 
of a charged D*, including contributions from B — > D**ir as well as non-resonant B — > D*inr 
decays; combinatoric B background, defined as all remaining BB background events; continuum 
e+e~ — ► qq, where q represents a, u, d, s, or c quark. 



4.3 Event Selection 

To suppress the continuum background, we select events in which the ratio of the 2nd to the 0th 
Fox- Wolfram moment [8], computed using all charged particles and EMC clusters not matched to 
tracks, is smaller than 0.40. Hard pion candidates are required to be reconstructed with at least 
twelve DCH hits. Kaons and leptons are rejected based on information from the IFR, DIRC, energy 
loss in the SVT and DCH, or the ratio of the candidate's EMC energy deposition to its momentum 
(E/p). We define the D* helicity angle 9d* to be the angle between the flight directions of the 
D and the B in the D* rest frame, calculated with = 0. Taking advantage of the longitudinal 
polarization in signal events, we suppress background by requiring |cos#£>*| to be larger than 
0.4. All candidates are required to be in the range 1.81 < m miss < 1.88 GeV/c 2 . When multiple 
candidates are found in the same event, only the one with the n?. m i SS value closest to M D o is used. 



4.4 Fisher Discriminant 

To further discriminate against continuum events, we combine fifteen event shape variables into 
a Fisher discriminant [9] F. Discrimination is provided due to the fact that qq events tend to 
be jet-like, whereas BB events have a more spherical energy distribution. Rather than applying 
requirements to the variable F, we maximize efficiency by using it in the fits described below. 
The fifteen variables are calculated using two sets of particles. Set 1 includes all tracks and EMC 
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clusters, excluding the hard and soft pion candidates; Set 2 is composed of Set 1, excluding all tracks 
and clusters with CM momentum within 1.25 rad of the CM momentum of the D, calculated with 
= 0. The variables, all calculated in the CM frame, are 1) the scalar sum of the momenta of all 
Set 1 tracks and EMC clusters in nine 20° angular bins centered about the hard pion direction; 
2) the value of the sphericity, computed with Set 1; 3) the angle with respect to the hard pion of 
the sphericity axis, computed with Set 2; 4) the direction of the particle of highest energy in Set 2 
with respect to the hard pion; 5) the absolute value of the vector sum of the momenta of all the 
particles in Set 2 ; 6) the momentum ph of the hard pion and its polar angle. 

4.5 Decay Time Measurement and Flavor Tagging 

We define z rcc to be the decay position along the beam axis of the partially reconstructed B 
candidate. To find z rec , we fit the hard pion track with a beam spot constraint in the plane 
perpendicular to the beams, the (x, y) plane. The actual vertical beam spot size is approximately 
5 ^m, but the constraint is taken to be 30 /im in the fit in order to account for the B flight distance 
in the (x, y) plane. The soft pion is not used in the fit, since it undergoes significant multiple 
scattering. 

The decay position z tag of the other B in the event (the tag B) along the beam axis is obtained 
from all other tracks in the event, excluding all tracks whose CM momentum is within 1 rad of the 
D CM momentum. The remaining tracks are fit with a beamspot constraint in the (x, y) plane. 
The track with the largest contribution to the \ 2 °f the vertex, if greater than 6, is removed from 
the vertex, and the fit is carried out again, until no track fails this requirement. 

We then calculate the decay distance Az = z rcc — z tag , and the decay-time difference At = 
Az/(7/3c). The machine boost parameter 7/3 is calculated from the beam energies, and its average 
value over the run period is 0.55. The vertex fits used to determine z rec and z tag also yield the Az 
error o^ z which is used to compute the event-by-event At error a At- 

The flavor of the tag B is determined from lepton and kaon candidates. The lepton CM momen- 
tum is required to be greater than 1.1 GeV/c in order to suppress "cascade" leptons originating in 
charm decays. We identify electron candidates using E/p, and the Cherenkov angle and number of 
photons detected in the DIRC. Muons are identified by the depth of penetration in the IFR. Kaons 
are identified using the ionization measured in the SVT and DCH, and the Cherenkov angle and 
number of photons detected in the DIRC. In either the lepton or kaon tagging category, if several 
tagging tracks are present, the track used for tagging is the one with the largest value of 9 tag-, the 
CM-frame angle between the track momentum and the D momentum. This is done in order to 
minimize the impact of tracks originating from the unreconstructed D. If there are both identified 
leptons and kaons in the event, the event is tagged using the lepton tracks only. 

We apply the following criteria in order to obtain good At resolution: the \ 2 probability of the 
Zrcc vertex fit must be greater than 0.001; at least two tracks must be used for the tag B vertex 
fit; 0"At is required to be less than 2 ps; and |At| is required to be less than 15 ps. To minimize 
the impact of tracks coming from the unreconstructed D, only tagging leptons (kaons) satifying 
cosOtag < 0.75 {cosOtag < 0.50) are retained. 

4.6 Probability Density Function 

The analysis is carried out with a series of unbinned maximum likelihood fits performed indepen- 
dently for the lepton-tagged and kaon-tagged events. The probability density function (PDF) is a 



11 



function of the missing mass m m i ss , the Fisher discriminant F, the decay time difference At, and 
its error a A t- 

The PDF for on-resonance data is a sum over the PDFs of the identified event types, 

P = $BB {fp+TT {fo-TT -Pd*tt + (1 - fD*w)PD*p) ^ 
+ (1 — fp+ir) /comb -Pcomb + (1 — /comb)-fpcakJ } + (1 — f Bb)^''M'> 

where is the PDF for events of type i, and fj are relative fractions of events, each limited to lie 
in the range [0, 1]. Each of the PDFs Pi is a product of the form 

Pi(m miss , F, At, a At , s tag , s mix ) = Mi(m miss ) Fi{F) T((At, a At, s tag , s m ix), (8) 

where the variables s tag and s m i x are determined by the flavor of the tag B and the charge of the 
hard pion 

Stag " \ -1, tag B = B° ' W 

J +1, unmixed , 1f> . 
Smix = { -1, mixed ' (10) 

and an event is labeled "unmixed" if the hard pion is a ir~(ir + ) and the tag B is tagged as a B°(B°) 
and "mixed" otherwise. The functions Mi, Ti, and T/ are described below. The parameters of Pi 
are different for each event type, except where indicated otherwise. 
The 

Timiss PDF for each event type i is the sum of a bifurcated Gaussian plus an ARGUS 
function: 

Mi{m miss ) = ff Gi(m miss ) + (1 - fi)Ai(m miss ), (11) 
where ff is the bifurcated Gaussian fraction. The functions Qi and Ai are 



G (r) oc {^[-{x-M t f/2al t ], x < M t 
A(x) oc x\jl- {x/Mff exp ^ (\ - (x/Mff^j 



(12) 



S(M t A -x), (13) 



where Mj is the peak of the bifurcated Gaussian, 0"^^ and am are its left and right widths, ti is 
the ARGUS exponent, Mf- is its end point, and the proportionality constants are such that each 
of these functions is normalized to unit area. 

The Fisher discriminant PDF Ti for each event type is a bifurcated Gaussian, as in Eq. 12. The 
parameter values of Td*^, Td* p , T^ cz ^, and T cou ^, are identical. 

The At-dependent part of the PDF for events of type i is a convolution of the form 

T((At, a A t, stag, Smix) = J dAt tTUC Ti(At truc , s tag , s mix ) Ki(At - At truc , a At ), (14) 

where Tj is the distribution of the true decay-time difference At true and IZi is a resolution function 
that accounts for detector resolution and effects such as systematic offsets in the measured positions 
of vertices. The resolution function for events of type i is the sum of three Gaussians: 

Ki(At - At truc , a At ) = f? g?(tr,<TM) + (1 " /I 1 " ft) 0T(t r , *At) + ff G°{t r , a At), (15) 
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where t r is the residual At — Ai true , and Q™, Qf, and Q° are the "narrow", "wide", and "outlier" 
Gaussians. The narrow and wide Gaussians have the form 



Gi(t r ,(TAt) = 



1 



2ir s-j a At 



■ exp 



/ (t r - b{a A t) 



4 VM? 



(16) 



where the index j takes the values j = n,w for the narrow and wide Gaussians, and bj and s\ are 
parameters determined by fits, as described in Sec. 5. The outlier Gaussian has the form 



Q° (t r , <7At) = 



1 



2vr s? 



exp 



2«) 2 



(17) 



where in all fits the values of b° and s° are fixed to ps and 8 ps, respectively, and are later varied 
to evaluate systematic errors. 

The PDF T D * n (At true ; Stag) Smix) for signal events corresponds to Eqs. 1 and 2 with Eq. 5 and 
additional parameters to account for imperfect flavor tagging. We define w^o (u)^ ) to be the 
mistag probability of signal events whose tag B was tagged as a B° (B°), when the tagging track is 
a daughter of the tag B. Then uj = (wgo + a - , j?o)/2 * s ^ ne avera § e niistag rate, and Aw = uj b o —w^o 
is the mistag rate difference. We further define a to be the probability that the tagging lepton or 
kaon is a daughter of the unreconstructed D produced in the B° — ► D^tt^ decay, and p to be the 
probability that this track results in a mixed flavor tag. With these definitions, the signal PDF is 
written as 

^~D*7r(Attruc) ^tag) ^mix) = 



Intone 

■ e r x 



{(l-a)[(l- Stag Aw) 

+«mix (1 - 2a;) cos(AmAttruc) 

-stag (1 - 2w) S ± sin(AmAt truc )] 

+a(l + *mix(l - 2p))} , (18) 

where the value ± in S ± is determined by the sign of the product s t agS m ix- The last term accounts 
for the tags due to daughters of the unreconstructed D. The parameters r, Aw, u, Am, and 
are determined from a fit to the data, as described below. 

The tag B may undergo a b — > ucd decay, and the kaon produced in the subsequent decay of the 
charmed meson may be used for tagging. This introduces additional terms, which are not present 
in Eq. 18. To take this effect into account, we use an alternative parameterization [10] for the kaon 
tags. In this parameterization [10], the coefficient of the sin(AmAt t rue) term in Eq. 18 changes, to 
give 



^~D*7r(Attrue) St a g, S m i x ) — . 



1 lAttruel 

e t x 

4r 



{(1-a) [(1- stag Aw) , 
+s m ix (1 - 2w) cos(AmAt truc ) 

- ((1 - 2 u) (s tag a + s mix c) + s t agSmix&(l ~ s tag Aw)) sin(AmAttruc)] 
+a(l + s mix (l - 2p))} , (19) 



where 



= 2r sin(2/3 + 7) cos S, 



(20) 
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b = 2r / sin(2/? + 7 )cos<5', (21) 

c = 2cos(2/3 + 7)(rsin(5- r'siaS'). (22) 

Here r' describes the effective ratio between the magnitudes of the b — > uca! and b — > cud amplitudes 
in the tag side decays, and <5' is the effective strong phase difference between these amplitudes. This 
parameterization neglects terms of order r 2 and r' 2 . 

We take the At PDF parameters of B — > D* T p ± events to be identical to those of the B° — > 
Z)* =F 7r ± events except for the CP- violating parameters S^, a, b, and c, which are set to and are 
later varied to evaluate systematic uncertainties. 

The Attruc PDF for the combinatoric and the peaking BB background have the same functional 
form as Eq. 18 but with independent values for the parameters. The parameterization of the At truc 
PDF for the peaking BB background has been determined from the Monte Carlo sample. 

The At truc PDF for the continuum background has the functional form 

i At truc 

T,5 = /* 5 <5(At true )7j 5 + (l-4)-— e ^ T£, (23) 

11 

where 

Tgg = 1 - Stag A^ + S mix (1-2^), 

T T m = 1 - stag A^ + s mix (1-2^). (24) 

5 ANALYSIS PROCEDURE 

The analysis takes place in four steps, each involving maximum likelihood fits, carried out simul- 
taneously on the on- and off-resonance data samples: 

1. Kinematic- variable fit: The parameters of -Mj(m m i ss ) and the value of /d*tt and in Eq. 7 
are obtained from the Monte Carlo simulation, conducted with the branching fractions from 
Ref. [7]. Using these parameter values, we fit the data using the PDF in Eq. 7, but with Eq. 8 
replaced by 

"i v '"'miss > 

F) = Mi(m miss ) Fi{F). (25) 

The parameters determined in this fit are f BB , f p +n, and / com b in Eq. 7, the parameters of 
-A4 g g(m m i ss ), and those of Fi(F) for both continuum and BB events. 

2. a and p fit: The kinematic-variable fit is repeated to determine the number of signal events 
above and below the cut on cos 6t ag (see section 4.5). These values are then used to compute 
the values of a and p in the At PDF (Eq. 18). This is done using values for the efficiencies 
of the cut on cos 9 tag determined from the Monte Carlo simulation. 

3. Sideband fit: We fit events in the m m i ss sideband 1.81 < m m i ss < 1.84 GeV/c 2 to obtain the 
parameters of the combinatoric BB PDF T^ omb (Attruej s tag> s m i x , <J/\t)- The PDF in Eq. 7 
is used in this fit, with f p + n = and / CO mb = 1, to account for the fact that the sideband 
is populated only by continuum and combinatoric BB events. The value of f BB - and the 
parameters of the continuum PDF T^(At true , s tag , s m i x , 0At) in the sideband are also floating 
in this fit. 
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4. Signal-region fit: Using the parameter values obtained in the previous steps, we fit the data 
in the signal region 1.845 < m m i ss < 1.880 GeV/c 2 . This fit determines all the floating pa- 
rameters of the signal PDF T^* 7r (At tr ue ) -Stag, s m i x , CAt), and the parameters of the continuum 
PDF T^(At t rue,stag,Smix,0-At) except for b° Dt7n 5^*^, 6°^, and sj^ of Eq. 17. 

In steps 3 and 4 we also fit for a possible difference between the B° and B° tagging efficiencies, 
which may be different for each event type. In order to minimize the possibility of experimenter 
bias, step 4 of the analysis is carried out in a "blind" manner, such that the values of S^, n are 
hidden from the analysts until all the systematic errors have been evaluated. 

The validity of the analysis procedure has been verified using the Monte Carlo simulation in 
two ways. First, the use of identical At PDFs for B — > D*^ p^ and B° — > D*^tt^ events (except 
for the CP violating parameters), as well as for the combinatoric BB background in the sideband 
and in the signal region, is validated by comparing the At distributions for these event types in 
Monte Carlo by means of Kolmogorov-Smirnov tests. The results of fits to the At distributions are 
also compared. In all cases, good agreement is observed. Second, the entire analysis procedure is 
carried out on a Monte Carlo sample containing four times the number of events observed in the 
data. The values of the parameters obtained in these Monte Carlo fits, most importantly, the CP 
parameters, are consistent with the input parameters to within the statistical uncertainties. In the 
case of the fit to the lepton-tagged events, a bias due to the assumption that events tagged with 
direct and cascade leptons are described by the same resolution function is studied using the full 
Monte Carlo simulation and a fast Monte Carlo technique. This bias is found to be +0.012 for 
Sp» n , and a corresponding correction is applied to the results presented in this paper. 

6 SYSTEMATIC STUDIES 

The systematic uncertainties on the CP violation parameters (S' ± for events tagged with a lepton 
candidate and a, b, and c for events tagged with a kaon candidate) are summarized in Table 1 and 
described here: 

(1) The statistical errors obtained in the kinematic- variable fit are propagated to the signal- 
region At fit. This is done by varying the parameters determined by the kinematic- variable 
fit, taking into account their correlated errors, repeating the signal-region At fit with the new 
parameters, and taking the resulting variation in the CP violation parameters as a systematic 
uncertainty. 

(2) The same procedure is applied for the statistical errors of the parameters determined in the 
sideband fit. 

(3) The parameters of the outlier Gaussian for the signal At PDF that are fixed in the signal- 
region fit are varied: s° is varied between 6 and 10 ps, and b° between —2 and 2 ps. The 
parameters a and p are varied by the statistical uncertainties determined in their fit. 

(4) We vary /d*tt according to the uncertainties in the branching fractions B(B° — > Z)*~7r + ) and 
B(B° — > D*~ p + ) [7], and repeat the signal region At fit. 

(5) In the signal-region At fit we set the values of Sp, to 0. To obtain the systematic error due 
to this, we vary S^, between —0.04 and +0.04. The same is done for S^ eak and <S^ mb and 
the resulting variations in the CP parameters are combined linearly to obtain a conservative 
uncertainty. 
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(6) The values of the neutral B meson lifetime r and the mixing frequency Am are left free in 
the fit. The fit is repeated setting Am to the published value [7]. It is also repeated with r 
set to the value obtained by fitting the Monte Carlo sample, which is lower than the value in 
Ref. [7] due to tracks originating from the unreconstructed D [4]. The resulting variations in 
the CP parameters are combined in quadrature to obtain the systematic uncertainty. 

(7) The parameters of the peaking background At PDF, nominally taken from Monte Carlo, 
are varied: we set their values to the values of the corresponding combinatoric background 
parameters, and alternatively fit the parameters w pea k) Am pea k, and T pea k from the data. The 
largest variation is used as a systematic error. We also fit signal Monte Carlo with and without 
BB backgrounds, taking the difference as a systematic error. An additional error accounts 
for a possible difference between the combinatoric BB parameters in the signal region and the 
sideband, evaluated by fitting Monte Carlo events with the parameters of T^ omb determined 
from BB Monte Carlo events in the signal-region or sideband. 

(8) The uncertainty due to the beam spot constraint is evaluated by performing the fit on a 
sample where the beam spot y position has been shifted by 20 fim. 

(9) The z length scale of the detector is determined with an uncertainty of 0.4% from the recon- 
struction of secondary interactions with a beam pipe of known length [11]. The uncertainty 
in the CP parameters due to the detector z scale is evaluated performing the fit on a sample 
where the z scale has been varied by 0.4%. 

(10) The effect of imperfect detector alignment on a time dependent CP violation measurement is 
estimated in a similar study using fully reconstructed decays and is assumed to be the same 
for this analysis. 

(11) The bias due to cascade leptons is evaluated by varying the parameters of the corresponding 
PDF. The fraction of cascade leptons is varied by 3% and the mistag rate by 2%. 

(12) The statistical error of the fit to the signal Monte Carlo sample is added to the systematic 
error to account for possible bias of the analysis procedure. 

7 RESULTS OF THE FITS 

The kinematic- variable fit yields 6406 ± 129 signal events for the lepton-tagged sample and 25157 ± 
323 signal events for the kaon-tagged sample. The results of the fit are shown in Figs. 3 and 4. The 
background in the lepton (kaon) sample is mainly due to BB (continuum) events. The kinematic- 
variable fit is then repeated on subsamples of the data, separated according to the flavor tag and 
the charge of the final state particles (Table 2). 

The fit to the signal region (Figs. 5 and 6), described in step 4 of Sec. 5, determines Am, r, u, 
Auj, and the CP parameters of the signal PDF TD* w (At tTue , Stag, Smix)- These parameters are S ± 
for lepton tags and a, b, c for kaon tags. Five parameters of the signal resolution function are also 
determined by the fit, as are eight continuum parameters: four parameters for the At tr ue PDF and 
four parameters for the resolution function. The results of the fits on the data are 

S+ = -0.078 ± 0.052 {stat.)± 0.020 (syst.), 

S- = -0.070 ± 0.052 (stat.) ± 0.021 (syst.) (26) 
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Table 1: The systematic uncertainties on the CP violation parameters: for events tagged with 
a lepton candidate, and a, b and c for events tagged with a kaon candidate. 



Source 


S error 


S + error 


a error 


b error 


c error 


(1) Kinematic- variable fit statistics 


0.0005 


0.0007 


0.0009 


0.0004 


0.0004 


(2) Sideband fit statistics 


0.0007 


0.0015 


0.0004 


0.0003 


0.0005 


(3) Variation of fixed At PDF parameters 


< io- 4 


< io- 4 


< 10~ 4 


< io- 4 


< io- 4 


(4) Uncertainty in branching fractions 


0.0028 


0.0028 


0.0017 


0.0002 


0.0033 


(5) Uncertainty in bkg CP parameters 


0.0096 


0.0096 


0.0128 


0.0073 


0.0127 


(6) Variation of r and Am 


0.0012 


0.0040 


0.0052 


0.0017 


0.0008 


(7) Effect of BB background 


0.0008 


0.0068 


0.0045 


0.0042 


0.0054 


(8) Beam spot 


0.0017 


0.0012 


0.0007 


0.0013 


0.0006 


(9) z scale 


0.0004 


0.0004 


0.0002 


< io- 4 


0.0003 


(10) Detector alignment 


0.0100 


0.0100 


0.0100 


0.0056 


0.0100 


(11) Cascade lepton bias 


0.0052 


0.0052 








(12) MC statistics 


0.0128 


0.0128 


0.0080 


0.0040 


0.0090 


Total 


0.020 


0.021 


0.019 


0.011 


0.020 



Table 2: Result of the kinematic-variable fit: number of signal events for the entire data sample 
and for samples separated according to the flavour tag and the charge of the final state particles. 





Lepton Tag 


Kaon Tag 


All 


6406 ± 129 


25157 ± 323 


B° tag 


3217 ±84 


12821 ± 232 


B° tag 


3179 ± 88 


12343 ± 224 


D*+7T- 


3136 ±85 


12299 ± 227 


D*~7T + 


3269 ± 84 


12830 ± 230 
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1.81 1.82 1.83 1.84 1.85 1.86 1.87 1.88 





-1.5 -1 -0.5 0.5 1 1.5 

Fisher Discriminant 




-1.5 -1 -0.5 0.5 1 1.5 



Fisher Discriminant 



Figure 3: Result of the kinematic- variable fit for lepton-tagged events. The (a) m m i ss and (b) Fisher 
discriminant F distributions for the on-resonance data are shown by the points with error bars. 
In the m. m i ss plot, the overlayed curves show, from bottom to top, the cumulative contributions of 
continuum, peaking BB, combinatoric BB, B° — > D* + p~ and B° — > D* + tt~ events. In the F plot, 
the PDFs for BB and continuum events are overlaid. Plots (c) and (d) show the same distributions 
for off-resonance data. 
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m miB ,(GeV/c) 



0.5 0.5 1 1.5 

Fisher Discriminant 




,(GeV/c) 



0.5 0.5 1 1.5 

Fisher Discriminant 



Figure 4: Result of the kinematic-variable fit for kaon-tagged events. The (a) m m i SS and (b) Fisher 
discriminant F distributions for the on-resonance data are shown by the points with error bars. 
In the m m i ss plot, the overlayed curves show, from bottom to top, the cumulative contributions of 
continuum, peaking BB, combinatoric BB, B — > D* + p~ and B — > D* + ir~~ events. In the F plot, 
the PDFs for BB and continuum events are overlaid. Plots (c) and (d) show the same distributions 
for off-resonance data. 



19 




A t (ps) A t (ps) 



Figure 5: Result of the At fit to the lepton-tagged data in the signal region. The points with error 
bars show the At distributions for (a) B°-tag unmixed, (b) -B°-tag mixed, (c) i?°-tag unmixed, (d) 
and i?°-tag mixed events. The curves show the PDF with the parameters obtained from the fit. 
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Figure 6: Result of the At fit to the kaon-tagged data in the signal region. The points with error 
bars show the At distributions for (a) B°-tag unmixed, (b) B°-tag mixed, (c) i?°-tag unmixed, (d) 
and i?°-tag mixed events. The curves show the PDF with the parameters obtained from the fit. 
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for the lepton-tagged events, and 

a = -0.054 ± 0.032 (stat.) ±0.019 (syst.), 
b = -0.009 ±0.019 (stat.) ±0.011 (syst.), 

c = ±0.005 ± 0.031 (stat.) ± 0.020 (syst.) (27) 

for the kaon-tagged events. We note the good agreement between a and (S + ± S~)/2. The largest 
correlation of S + (S~) with any linear combination of the other parameters floated in the fit is 0.24 
(0.26) and the correlation between S + and S~ is —0.057. 

We define two time-dependent CP-violating asymmetries from the numbers of events observed 
at time t with specific combinations of flavor tag and reconstructed final state: 

JV(tag B°, D* ± 7r =F )(t) - JV(tag B°, D^ir^jt) 
iV(tag B ,D*±7r^)(t) + N(tag B°, D*±^j(J~) ' 

iV(tag B° ± B°, D*-7r+)(t) - iV(tag B° ± B°, D* + TT~)(t) 
N(tag B° ± B°, D*-TT+)(t) ± JV(tag B° ± B°, D*+ir~)(t)' 

In the absence of background and experimental effects, A r Jp = — 2r sin(2/3 ± 7) cos 5 sin(AmAt) 
and Aqp = 2r'sin(2/3 ± 7) cos 5' sin (Am At). The asymmetry plots obtained for the data in the 
signal region are shown in Fig. 7 for lepton tags and in Fig. 8 for kaon tags. As expected, no 
time-dependent asymmetry is visible for A^p in the lepton case. 



Arec 
JK CP 



A^p 



8 PHYSICS RESULTS 

Combining a and (S + ± S~~)/2 from Eqs. 26 and 27, accounting for correlated errors, we obtain 

2r sin(2/3 ± 7) cos 5 = -0.063 ± 0.024 (stat.) ± 0.017 (syst.). (29) 

This measurement deviates from by 2.1 standard deviations, and is the main result of this analysis. 
From the difference (S + — S~)/2, we obtain 

2r cos(2/3 ± 7) sin S = -0.004 ± 0.037 (stat.) ± 0.020 (syst.). (30) 

We use two different methods for extracting constraints on | sin(2/3 ± 7)) from our results. 
We emphasize that the two methods make use of different additional information and different 
assumptions, and are therefore not directly comparable. These constraints are interpretations of 
our experimental results. Each of the methods involves defining and minimizing a x 2 function of 
sin(2/3±7) and other parameters. The x 2 functions are symmetric under the exchange sin(2/3±7) — > 
— sin(2/3 ± 7). Due to the large uncertainties and the fact that the minimum value of the x 2 may 
occur at the boundary of the physical region (| sin(2/3 ± 7)) = 1), the errors naively obtained from 
the variation of the x 2 functions are not relevant. In order to give a probabilistic interpretation to 
the results, we apply the Feldman-Cousins method [12] to set limits on | sin(2/3 ±7)]. 

In the first method we make no assumption regarding the value of r and use no additional 
experimental information about r. In this method, for different values of r we minimize the function 

3 

X 2 = E AxjVfSAxk, (31) 
j,k=i 
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Figure 7: Raw asymmetries Aqp and A^p 
(Eq. 28) for lepton tags as a function of At. The 
curves show the projection of the fitted PDF. 



Figure 8: Raw asymmetries A r <jp and A^p 
(Eq. 28) for kaon tags as a function of At. The 
curves show the projection of the fitted PDF. 
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where Axj is the difference between the result of our measurement and the theoretical expression 
for S + (j = 1), S~ (j = 2), and a (J = 3), and V is the measurement error matrix, which is 
almost diagonal. The parameters determined by this fit are sin(2/3 + 7), which is limited to lie in 
the range [—1, 1], and 5. The measurements of b and c are not used in the fit, since they depend 
on the unknown values of r' and 5'. We generate many parameterized MC experiments with the 
same sensitivity as reported here for different values of sin(2/3 + 7) and r. The fraction of these 
experiments in which x 2 (sin(2/?+7)) — Xmin * s smaller than in the data is computed and interpreted 
as the confidence level (CL) of the lower limit on | sin(2/5 + 7)|. The 90% and 95% CL limits as a 
function of r are shown in Fig. 9. The x 2 fit determines | sin(2/3 + 7)] up to the twofold ambiguity 
I sin(2/3 + 7)| <-> I cos<5|. The limits shown in Fig. 9 are always the more conservative of the two 
possibilities. 
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Figure 9: The 90% and 95% CL lower limits on | sin(2/3+7)| as a function of r, using no experimental 
information on r. 

In the second method we assume that r may be estimated from 



r = w iBtp^D-jTlfr (32) 
6 V B(B« -. D--W+) f D; 1 ' 

where 9q is the Cabibbo angle. Using the branching fractions B(B° — > D*~n + ) = (0.276 ± 
0.021)% [7], B(B° -> D*fir-) = (1.9+^3 ± 0.5) x lO" 5 [13] and the ratio of decay constants 

j& = 1.10 ±0.02 [14] yields 

r = 0.017^7- (33) 

An additional non-Gaussian 30% relative error is associated with the theoretical assumptions in- 
volved in obtaining this value. To carry out this method, we minimize 

3 

X 2 = £ Ax^Azfc + AV), (34) 
j,k=i 
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Figure 10: x 2 — Xmin as a function of sin(2/3 + 7), assuming r 
non-Gaussian theoretical uncertainty. 
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Figure 11: Confidence level (CL) as a function of the lower limit on | sin(2/3 + 7)], computed as 
the fraction of parameterized MC experiments where x 2 (sin(2/3 + 7)) — x m in 1S l ar ger than in the 
data, with the assumptions of Fig. 10. The vertical lines correspond to 90% and 95% confidence 
level and lead to the bound | sin(2/3 + 7 )| > 0.75 (0.62) at the 90% (95%) CL. 
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where the function 



A 2 (r) 



' ( r-r y 
V 0.005/ 




(r - r )/r > 0.3, 
k - r o|/?*o < 0.3, 
limmY (r-r )/r <-0.3 



(35) 



accounts for the 30% theoretical error and the Gaussian experimental error around the central 
value r = ro, Eq. 34. In addition to sin(2/3 + 7) and 6, the parameter r is also determined by this 
fit. The minimum of % 2 occurs at | sin(2/3 + 7)! = 1.0, 5 = 0, and r = 0.024, and at a value of 
Xmin = 0-54 for one degree of freedom. The value of x 2 — Xmin as a function of sin(2/3 + 7) and the 
resulting Feldman-Cousins confidence level curve are shown in Figs. 10 and 11. This method yields 
the limits | sin(2/3 + 7 )| >0.88 at 68% CL, | sin(2/3 + 7 )| >0.75 at 90% CL, and | sin(2/3 + 7)| >0.62 
at 95% CL. 



9 SUMMARY 

We present preliminary results of a study of time-dependent CP asymmetries in the B° — > I?* T 7r ± 
decay channels using the partial reconstruction method. The time-dependent CP asymmetry that 
we measure, 

2r sin(2/3 + 7) cos S = -0.063 ± 0.024 (stat.) ± 0.017 (syst.), (36) 

is different from by 2.1 standard deviations. This asymmetry does not depend on assumptions 
regarding r, the ratio of the magnitudes of the b — ► ucd and b — > cud amplitudes contributing to 
this decay. We present model-independent bounds on | sin(2/3 + 7)] as a function of r. With some 
assumptions regarding r, our results can be interpreted as a limit on the combination of CKM 
angles 2/3 + 7, | sin(2/3 + 7 )| > 0.75 (0.62) at the 90% (95%) CL. 
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